**************************************************************************************************
*** This do file creates the replication results for 
*** Eric Neumayer and Thomas Pluemper
*** Does ‘Data Fudging’ Explain the Autocratic Advantage? Evidence from the Gap between Official Covid-19 Mortality and Excess Mortality
*** Published in: Social Science & Medicine – Population Health (forthcoming) https://doi.org/10.1016/j.ssmph.2022.101247  
							
**************************************************************************************************
**************************************************************************************************
/* Note: You have to change "local DIR" to the directory you copy the original stata files contained 	*/
/* in the zip file and then run the do file. 
/* You must have Stata version 13 or higher installed													*/
**************************************************************************************************
*/

drop _all
clear matrix
clear mata


capture net install outreg2, from(http://fmwww.bc.edu/RePEc/bocode/o) /* checks whether outreg2 is installed 		*/



***********************************************************************************
local DIR = "D:\Research\Health for all database\Europe"  /*change relative path to the directory where the files are located */
cd "`DIR'"
***********************************************************************************




*** table 2

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)


* Covid-19, reported only
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, replace excel


* Covid-19, all
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel

* Excess mortality, reported only
reg excess_deaths_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel
 

* Excess mortality, all
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


* Gap, reported only
reg DV_excess_ab_C19_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel

* Gap, all
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


*** Robustness tests

** Freedom House measure of democracy

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  e_fh_pr_2020 i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)

* Covid-19, reported only
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_fh_pr_2020  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, replace excel


* Covid-19, all
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_fh_pr_2020  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel

* Excess mortality, reported only
reg excess_deaths_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_fh_pr_2020  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel
 

* Excess mortality, all
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  e_fh_pr_2020 i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


* Gap, reported only
reg DV_excess_ab_C19_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_fh_pr_2020   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel

* Gap, all
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_fh_pr_2020  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel

** EIU measure of democracy

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  e_fh_pr_2020 i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)

* Covid-19, reported only
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 democracy_eiu  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, replace excel


* Covid-19, all
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 democracy_eiu  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel

* Excess mortality, reported only
reg excess_deaths_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 democracy_eiu  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel
 

* Excess mortality, all
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  democracy_eiu i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


* Gap, reported only
reg DV_excess_ab_C19_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 democracy_eiu   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel

* Gap, all
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 democracy_eiu  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel



** Polity measure of democracy

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  e_fh_pr_2020 i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)

* Covid-19, reported only
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_polity2_2018  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, replace excel


* Covid-19, all
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_polity2_2018  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel

* Excess mortality, reported only
reg excess_deaths_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_polity2_2018  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel
 

* Excess mortality, all
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  e_polity2_2018 i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


* Gap, reported only
reg DV_excess_ab_C19_per_100k  new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_polity2_2018   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
outreg2 using table1, append excel

* Gap, all
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 e_polity2_2018  i.northern_hemisphere##i.periodvar, cluster(country)
outreg2 using table1, append excel


** Higher order polynomials of control variables

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020  v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)


* Covid-19, reported only
reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)

* Covid-19, all
reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)

* Excess mortality, reported only
reg excess_deaths_per_100k  c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)

* Excess mortality, all
reg excess_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020  v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)

* Gap, reported only
reg DV_excess_ab_C19_per_100k  c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 v2x_libdem   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)

* Gap, all
reg DV_excess_ab_C19_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)


** Higher order polynomials of all variables, including democracy

use "Article for SS&M PH (C19)", clear

gen lstringency_index=l.stringency_index

quietly reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020  c.v2x_libdem##c.v2x_libdem##c.v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)
keep if e(sample)


* Covid-19, reported only
reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 c.v2x_libdem##c.v2x_libdem##c.v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))

* Covid-19, all
reg covid_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 c.v2x_libdem##c.v2x_libdem##c.v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))

* Excess mortality, reported only
reg excess_deaths_per_100k  c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 c.v2x_libdem##c.v2x_libdem##c.v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))

* Excess mortality, all
reg excess_deaths_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020  c.v2x_libdem##c.v2x_libdem##c.v2x_libdem i.northern_hemisphere##i.periodvar, cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))

* Gap, reported only
reg DV_excess_ab_C19_per_100k  c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 c.v2x_libdem##c.v2x_libdem##c.v2x_libdem   i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))

* Gap, all
reg DV_excess_ab_C19_per_100k c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million##c.new_cases_smoothed_per_million c.median_age##c.median_age##c.median_age c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund##c.people_fully_vaccinated_per_hund  c.ln_gdp_per_capita##c.ln_gdp_per_capita##c.ln_gdp_per_capita   c.lstringency_index##c.lstringency_index##c.lstringency_index c.e_wbgi_gee_2020##c.e_wbgi_gee_2020##c.e_wbgi_gee_2020 c.v2x_libdem##c.v2x_libdem##c.v2x_libdem  i.northern_hemisphere##i.periodvar, cluster(country)
margins, dydx(v2x_libdem) at(v2x_libdem=(0 (0.15) 0.9))


** Yearly instead of monthly data
* Reported only
use "Article for SS&M PH (C19) yearly reported", clear

gen lstringency_index=stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem , cluster(country)
keep if e(sample)

* Covid-19, reported
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  , cluster(country)
outreg2 using table1, replace excel


* Excess mortality, reported
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem , cluster(country)
outreg2 using table1, append excel

* Gap, reported
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  , cluster(country)
outreg2 using table1, append excel

* All
use "Article for SS&M PH (C19) yearly all", clear

gen lstringency_index=stringency_index

quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem , cluster(country)
keep if e(sample)

* Covid-19, all
reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  , cluster(country)
outreg2 using table1, replace excel


* Excess mortality, all
reg excess_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020  v2x_libdem , cluster(country)
outreg2 using table1, append excel

* Gap, all
reg DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  , cluster(country)
outreg2 using table1, append excel



*** Appendix 1
use "Article for SS&M PH (C19)", clear
capture drop lstringency_index
gen lstringency_index=l.stringency_index

* Covid-19, full
quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
keep if e(sample)
collapse year, by(country)
drop year
list

use "Article for SS&M PH (C19)", clear
capture drop lstringency_index
gen lstringency_index=l.stringency_index

* Covid-19, full
quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Predicted" , cluster(country)
keep if e(sample)
collapse year, by(country)
drop year
list


*** Appendix 2
use "Article for SS&M PH (C19)", clear
capture drop lstringency_index
gen lstringency_index=l.stringency_index

* Covid-19, full
quietly reg covid_deaths_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita   lstringency_index e_wbgi_gee_2020 v2x_libdem  i.northern_hemisphere##i.periodvar if type=="Reported", cluster(country)
keep if e(sample)
sum covid_deaths_per_100k excess_deaths_per_100k DV_excess_ab_C19_per_100k new_cases_smoothed_per_million median_age people_fully_vaccinated_per_hund  ln_gdp_per_capita  lstringency_index e_wbgi_gee_2020 v2x_libdem
